Stabilizing a Bell state of two superconducting qubits by dissipation engineering 
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We propose a dissipation engineering scheme that prepares and protects a maximally entangled 
state of a pair of superconducting qubits. This is done by off-resonantly coupling the two qubits 
to a low-Q cavity mode playing the role of a dissipative reservoir. We engineer this coupling by 
applying six continuous- wave microwave drives with appropriate frequencies. The two qubits need 
not be identical. We show that our approach does not require any fine-tuning of the parameters and 
requires only that certain ratios between them be large. With currently achievable coherence times, 
simulations indicate that a Bell state can be maintained over arbitrary long times with fidelities 
above 94%. Such performance leads to a significant violation of Bell's inequality (CHSH correlation 
larger than 2.6) for arbitrary long times. 



Entanglement is a fundamental, yet counter-intuitive 
concept in quantum mechanics. Maximally entangled 
two-qubit states, often called Bell states, violate classical 
correlation properties [iHl] and are an essential build- 
ing block for quantum communication and quantum in- 
formation. Unfortunately, these states are also difficult 
to generate and sustain as interaction with the environ- 
ment typically leads to rapid loss of their unique quantum 
properties. Therefore, stabilizing a Bell state is a sought 
after goal. 

Quantum state stabilization can be achieved by an ac- 
tive feedback loop in which the system is measured, and 
conditioned on the measurement of an error, a gate re- 
stores the system to the desired state [SJ This ap- 
proach suffers from the latency of data acquisition and 
analysis that must take place during the lifetime of the 
quantum system, a particularly acute problem for rel- 
atively short-lived quantum systems such as supercon- 
ducting qubits. An alternative approach to stabiliza- 
tion is quantum dissipation (reservoir) engineering [7|. 
In this approach, autonomous feedback is achieved by 
drive-generated couplings between the quantum system 
and a dissipative reservoir |51fTS]. Thus, this type of au- 
tonomous feedback corrects errors without any need for 
real-time processing of the measured record. 

In the field of quantum Josephson circuits, autonomous 
stabilization of single qubit states has recently been im- 
plemented by using a low-Q cavity mode dissipa- 
tive reservoir |17l [T5| . While schemes for autonomous 
stabilization of entangled states exist in atomic physics 
contexts [SJ [TDl EH HI], the required symmetry in the 
physical system would appear to be an obstacle to im- 
plementing these with superconducting qubits. Indeed, 
while qubits in atomic physics experiments naturally pos- 
sess identical frequencies and coupling strengths to an 
engineered environment, such symmetry is difficult to 
achieve with superconducting artificial atoms. On the 
other hand, since these systems are engineered, they ben- 
efit from a large design versatility. In circuits, a large set 
of Hamiltonians can be designed for dissipation engineer- 



ing schemes. Moreover, coupling strengths can be tuned 
by simple modifications of circuit elements, and hence 
are not fixed by fundamental constants, as in the case of 
atomic physics systems. 

In this Letter, we propose to stabilize a maximally en- 
tangled state of a pair of superconducting qubits. Since 
we require no strong built-in symmetry in the physi- 
cal setup, the proposed method is readily implementable 
with superconducting qubits and should lead to break- 
through performance with currently achievable exper- 
imental parameters. The scheme exploits the strong 
dispersive interaction |19| between two superconduct- 
ing qubits and a single low-Q cavity mode and only 
requires continuous-wave (CW) microwave drives with 
well-chosen fixed frequencies, amplitudes and phases (see 
Fig. [TJ. Moreover, our scheme is robust against small 
variations of the latter control parameters, and requires 
only some basic calibration. By avoiding resonant in- 
teractions between the qubits and the cavity mode the 
qubits remain protected against the Purcell effect which 
would reduce their coherence times. 

Here, the low-Q cavity mode acts as an engineered 
reservoir which evacuates entropy from the qubits when 
a perturbation occurs: by driving the qubits and cav- 
ity with CW drives, we induce an autonomous feedback 
loop which corrects the state of the qubits every time it 
decays out of the desired Bell state. In order to show 
the robustness of our scheme and its ability to prepare 
and protect a state significantly violating Bell's inequal- 
ity, we have performed simulations with experimentally 
achievable parameters and realistic decoherence values. 

Consider two non-identical qubits coupled to a low-Q 
cavity mode (see Fig. [l|. The Hamiltonian of the three 
body system (two qubits and one cavity mode) in the 
Jaynes-Cummings approximation reads [20j : 

H = ^ + ^ + w^aTa, 
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FIG. 1: Possible setup for the proposed experiment: two su- 
perconducting qubits (Alice and Bob) interact with a mode 
of a 3D cavity. CW tones are applied on the input port of 
the cavity. These tones drive the qubits and the cavity and 
stabilize a Bell state of the two qubits. This stabilization is 
autonomous: no analysis of the output signal is needed. 



where cr^, erf are the PauH a ^ operators for the Al- 
ice and Bob qubit respectively, and a^^(a) is the cav- 
ity creation (annihilation) operator. Here, w^, ws and 
ijjc are the resonance frequencies for Alice, Bob and the 
cavity mode respectively, and g^, are the coupling 
constants of the cavity to Alice and Bob respectively. 
In the dispersive regime, \^a\ = |wa — wd ^ gA and 
I As I — \lob ~ Wc| S> 5s, and in the rotating frame for 
both qubits and the cavity, we obtain the following effec- 
tive Hamiltonian: 




FIG. 2: (a), (c) Spectra of qubit Alice and the cavity, respec- 
tively. The peaks are resolved since xa,Xb > k, I/T2 ■ 
(b) Energy level diagram of two qubits off-resonantly coupled 
to a cavity mode. Straight line arrows indicate couplings be- 
tween energy levels which are induced by CW drives. The 
wavy arrow indicates energy decay from the cavity mode. 
These drives couple the two qubits to the cavity, here used 
as a quantum reservoir, in such a way that the Bell state 
\4>-) = (Iffe) — |e3))/\/2 is the only stable state. 
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where xa = '^9/? I and xs = '^gs^/^B are the dis- 
persive coupling strengths between the qubits and the 
cavity mode. We neglect the direct qubit-qubit inter- 
action, which has no effect on our stabilization proto- 
col. The same effective Hamiltonian would be obtained 
by more realistic quantum circuit models [2T]. Here, 
we consider the coupling to be in the strong dispersive 
regime |19| . where xa and xb are assumed to be larger 
than the cavity and qubit line- widths xa,Xb > k ^ 
l/r^,l/T^. In this regime, the resonance frequency of 
each qubit depends on the number of photons in the cav- 
ity, and conversely the cavity frequency depends on the 
states of the qubits (see Fig. |2] [a,c]). We can, hence, 
drive the qubits (cavity mode) conditionally on the state 
of the cavity mode (qubits). The stabilization drives 
lead to a measurement induced dephasing at a rate 
which increases with the difference of the dispersive shifts 
\xa — Xb\- Thus, the only symmetry requirement for our 
dissipation engineering scheme is that xa, Xb are suffi- 
ciently close to each other to ensure Km ^ since k is 
the time scale associated with the autonomous feedback 
loop. Using |[22j, we show that this condition is satisfied 
when 



\XA - XbI < XaXb/kVti 



(2) 



Note that such symmetry is easily achieved with super- 
conducting qubits, by tuning the frequency of one of 
the qubits so that 2gA'^ / ^a and 2gB^/^B become close 
enough. 



We denote as wf^, wf^, and w^^ the cavity fre- 
quencies when the qubits are respectively in the states 
1 35) 1 1 ^e) , leg) and |ee). We have used the notation 
Iffs) — Ifi')'^' |e), where the first (second) element refers to 
the state of the Alice (Bob) qubit. Similarly we denote 
u!^ {ujg) the frequency of Alice (Bob) when there are n 
photons in the cavity. 

We now detail, step by step, the dissipation engineering 
procedure by explaining the task performed by each CW 
drive (see Fig. |2]). 

Parity selection: Applying two CW drives on the 
resonator at frequencies cjf^ and w^'^ with equal ampli- 
tudes Ec = f V^i we generate a coherent state with mean 
photon number n in the cavity when the qubits are in 
\gg) or |ee). Hence, an initial state \gg,0) (|ee, 0)) will 



,- XA+XE 



converge to a state gg, ae 

where [ae*"^) is a coherent state with 
other hand, states of the form 



ee, ae 



: XA+XE 



left invariant 
Xa,Xb > e, 



Indeed, 



-eg 



\eg,0) 



- n. On the 

;ge l^e, 0) are 



we assume \lo'^^^'^^ — w^'^'^fj — 
and hence the drives are off- 
resonant with the cavity when the qubits are in an odd 
parity state, leaving the cavity in vacuum. Note that 
for finite ratios xa.b/^c, the cavity will populate slightly 
when the qubits are in an odd parity state, inducing a 
dephasing which will be much smaller than the cor- 
rection rate n as long as condition ^ is satisfied. At this 
stage, if we measure the photons leaking out of the cavity, 
we can determine if the state of the qubits is \gg) , |ee) or 
in the odd parity manifold spanjlgie) , leg)}. 

Bell state selection: We now apply two Rabi drives 



3 



separating the two Bell states \4>-) = {\ge) — \eg))/y/2 
and \4>+) = (l^e) + \eg))/y/2 in the odd parity mani- 
fold. Since the cavity mode remains in the vacuum state 
within this manifold, we apply two qubit drives at fre- 
quencies u}\ and Wg, and furthermore the corresponding 
Rabi amplitudes are chosen to be the same and given by 
n". Assuming with no loss of generality that these drives 



are around the X axis, this adds the term r2°(c 



Oto 



the Hamiltonian in the rotating frame of the qubits. The 
latter term leaves invariant the state | </>_), and couples 
the state \4>+) to [\gg) + \ee)) / ^/2. These two qubit drives, 
added to the two cavity drives, pump the population from 

|(/)-(-,0) to a mixture of the states gg,ae^ 

XA + Xl 



and 



ee,ae 2 -y^ while keeping the state |(/)_,0) invari- 
ant. Hence, by looking at the photons leaking out of the 
cavity, we can determine if the qubits are in the desired 
target state {(j)-), or if they have been projected to \gg) 
or |ee). 

Irreversible pumping to the target Bell state: 

To complete the dissipation engineering scheme which 
stabilizes \4>-), we need to pump the population from 



gg,ae- 



and 



ee, ae 



: XA+XE 



back to 



To 



this end, we apply two more drives at frequencies 



; XA+XB 



and 



+ n 



XA+XE 



with the same Rabi fre- 



quency fi". Adjusting the phase difference of these 
drives appropriately, we obtain the following term in 
the Hamiltonian (in the rotating frame of the qubits): 
il" |^g-2«''-*2*^"° *("■+ ^ cr+) +c.c.^ (c.c. stands for com- 
plex conjugate). For integer n Ri n, these drives couple 
the states \gg,n) and \ee,n) to |0_,n). The latter state 
is unaffected by the cavity drives and so decays down to 
\(j}-,0), as indicated by Fig. [2] 

In summary, turning on the above six CW drives simul- 
taneously pumps the population from any of the states 
\gg), |ee) and |0+) to the state 10-). This pumping hap- 
pens at a rate of order k: the states \gg, 0) and |ee, 0) con- 



verge to the states ^(7, ae'"'* 2''^ *^ and ee, ae 

at rate k; the state \(j)_,n) decays towards |0_,O) at rate 
k; and we can choose the oscillation rates il'^ and ft" at 
will. Provided that this pumping rate of order k is much 
larger than the decoherence rates of the qubits that de- 
stroy entanglement, the dissipation engineering efficiently 
prepares and protects the desired Bell state. Note, fur- 
thermore, that any imperfection in the experimental pa- 
rameters, such as non-ideal tuning of the drives and a 
difference between XA and xb can also lead to loss of 
coherence for both qubits: however, as long as the in- 
duced dephasing rate is smaller than the pumping rate, 
the protocol remains efficient. 

We present simulations of the Lindblad master equa- 
tion with practically achievable parameters to address 
all these imperfections and demonstrate our scheme is 
robust to them. The system parameters we take for 
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FIG. 3: Simulation of the fidelity with respect to the tar- 
get state \<j>-} as a function of time. The CHSH correla- 
tion associated to each fidelity can be read on the right axis. 
The system parameters are Tj^j^ — 50 ^s, K,/2Tr — 2 MHz, 
Xa/2tv = 10 MHz, Xb/2tv = 9.5 MHz. The cavity is driven 
with ft — 4, the qubit drive strengths are il" = | and Q," = k 
and the initial state is the ground state \gg,0). Notice that 
realistic decoherence mechanisms and imperfect parameter 
tuning are inciuded in the simulation. Despite these imper- 
fections, the desired Bell state is stabilized with a fidelity of 
about 94% resulting in a constant violation of Bell's inequal- 
ity (stable CHSH correlation of 2.64 close to the maximal 
violation limit 2\/2 ^ 2.83). 



the simulation are 



A,B 



50 ^s, k/2t: ^ 2 MHz, 



Xa/Stt = 10 MHz, xs/27r = 9.5 MHz. For transmon 
qubits, a Purcell filter j23] might be necessary to obtain 
KT^f » 1. 

We simulate the Lindblad master equation 



dt 



Pit) 



[H(t),p(t)]+«D[a]p(t) 




where 



H(t) = ixA- 



XB^)a^a 



r, iXA+ Xb ru , u 

2ecC0s( t){a + aJ) 



2 

, XA+XE 



c.c. 



^A,B 



with cr+ = ((T_)^^ = |e) (g|. Furthermore, 1/T^ 

— 1/2T^'^ are the pure dephasing rates and the 
Lindblad super-operator is defined for any observable O 



D[0]p = OpO^ - -O^Op - -pO^O 
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Fidelity to Bell state 




FIG. 4: Fidelity (color) of the two qubit steady state with 
respect to the pure target state {(j}-), as a function of the two 
drive strengths. The size of the 90% fidelity contour demon- 
strates the robustness of the protocol. The parameter domain 
where a Bell violation would occur is even larger (75% con- 
tour). Large n contradicts the constraint ||2|, thus leading to 
a higher induced dephasing between the two qubits. On the 
other hand, for small n, we are constrained to a small range 
of in order to satisfy nxA,B ^ which is needed to per- 
form conditional qubit operations. Taking the parameters of 
the previous simulations, the optimal compromise is around 
3 < n < 5 and 0.5 < / k < 1, leading to a fidelity greater 
than 93%. 



Choosing optimal drive amplitudes, we present the evo- 
lution of our system in Fig. |3] In this figure we plot the 
fidelity F{t) of the two qubit state with respect to the 
pure target state \(j>-). The corresponding CHSH cor- 
relation B{t) can be read on the right hand side axis. 
More precisely, we have F{t) = Tr ((|0_) {<j>- \ (8) Ic)p{t)) 
and Bit) = Tr ((Ochsh ® Ic)p(t)), where 



O 



CHSH 



V2 



V2 



A'^y 



V2 



72 ' 



and Ic is the identity operator on the cavity. The state 
violates Bell's inequality when the CHSH correlation rises 
above 2 (dash-dotted green horizontal line) [U] . 

Fixing the experimental parameters, one can scan the 
drive amplitudes to find the maximal violation of the Bell 
inequality. This has been performed in Fig. |4] where we 
scan n and il". As can be seen, for a large range of ampli- 
tudes, the Bell inequality remains violated which shows 
the robustness of our protocol to these control parame- 
ters. A discussion on the optimal values of the amplitudes 
is provided in the caption of the figure. 

In conclusion, we have presented a scheme to stabilize a 
Bell state without the need for active feedback, but rather 
by engineering a coupling to a dissipative reservoir which 



performs an autonomous feedback loop. Our protocol is 
robust to system and control parameters, and does not 
require the two qubits to be identical. Numerical simu- 
lations taking into account all relevant experimental im- 
perfections predict a CHSH correlation above 2.6 (higher 
than 94% fidelity to the target Bell state) clearly violat- 
ing Bell's inequality for arbitrarily long times. This Bell 
state could be further purified by monitoring the cavity 
output and selecting the state only when no photons are 
leaking out from the cavity. Preliminary experimental 
results are in excellent agreement with this theory and 
will be published elsewhere. 
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